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Equations of motion are derived for a general spacecraft composed of rigid bodies con- 
nected via rotary (spherical or gimballed) joints in a tree topology. Several supporting 
concepts are developed in depth. Basis dyads aid in the transition from basis-free vector 
equations to component- wise equations. Joint partials allow abstraction of 1-DOF, 2-DOF, 
3-DOF gimballed and spherical rotational joints to a common notation. The basic building 
block consisting of an “inner” body and an “outer” body connected by a joint enables ef- 
ficient organization of arbitrary tree structures. Kane’s equation is recast in a form which 
facilitates systematic assembly of large systems of equations, and exposes a relationship 
of Kane’s equation to Newton and Euler’s equations which is obscured by the usual pre- 
sentation. The resulting system of dynamic equations is of minimum dimension, and is 
suitable for numerical solution by computer. Implementation is discussed, and illustrative 
simulation results are presented. 


Nomenclature 

Reference Frames and Points 

N A Newtonian (non-accelerating, non-rotating) frame 

B A Body frame 

Bj The jth Body frame 

B* Mass center of Bj 

Gk The kth Joint 

Bi Inner Body 

B a Outer Body 

Vector-Dyad Notation 
x A vector 

Xj The components of x in frame Bj 
B_ The Basis Dyad for frame B 
Partial Angular Velocity Dyad 
Y_ Partial Velocity Dyad 
to Angular velocity of B in N 
v Linear velocity of B* in N 

bk Orthonormal basis vectors for B, ( k = 1,2,3) 
bjk Orthonormal basis vectors for Bj, ( k = 1,2,3) 
hk Orthonormal basis vectors for N, ( k = 1, 2, 3) 

Fjk The path vector from Bj to joint Gk 

Pj The path vector from B* to B\ 

Matrix Notation 
U Identity Matrix 

0 Partial Angular Velocity Matrix 

V Partial Velocity Matrix 
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u Dynamic states (generalized speeds) 

x Kinematic states (generalized coordinates) 

w Components of uj, expressed in B 

v Components of F, expressed in N 

6 Gimbal angles 

a Joint angular rates 

T Joint partial matrix 

H Angular momentum 

m Mass 

I Central Moments of Inertia 

F Force 

T Torque 

a Acceleration 

a Angular Acceleration 

a T Remainder acceleration 

a r Remainder angular acceleration 

iC k Direction Cosine Matrix between frames Bj and Bk 

r* The skew-symmetric matrix such that r x x = f x x 

I. Introduction 

A broad range of spacecraft may be well modeled as a collection of rigid bodies connected in a tree 
topology by rotational joints. The complexity of appropriate dynamical models can quickly outstrip the 
limited examples available in textbooks. Powerful and general commercial dynamics software packages are 
available, but require an investment of money and training time, and present the challenge of integrating 
their dynamics capabilities with custom environmental and control models. Consequently this practitioner, 
like many before, undertook to develop custom simulation code to solve the range of problems that we face in 
our daily work. The work presented here is one practitioner’s journey from textbook to working simulation 
code. There is no “new” theory introduced, but we emphasize some concepts that we found key to successful 
application. The objective is to obtain a system of differential equations of motion amenable to solution 
by numerical integration. We talk about Kane’s equation, degrees of freedom and constraints, basis dyads, 
partial velocities and joint partials. We focus on the building block composed of a single joint connecting an 
“inner” and “outer” body, then use that to assemble trees of arbitrary complexity. We present a top-level 
algorithm of our simulation code, and present some illustrative simulation results. 


II. Kane’s Equation 

Kane’s equation is presented as a set of scalar equations: 1 

F r + Ff = 0, r = l,...,n 


(1) 


where n is the number of degrees of freedom of the system. For a system comprising a set of N b rigid bodies, 
the generalized active forces are given by 


N b N b 

F r = 5^ ft) +$>;.&) 

fc=l k= 1 


and the generalized inertia forces are given by 

N b N b 

f; = £ [af • (-as* - <3* x S k )] + [a? • (-m*ay] 

k = 1 


k = 1 

Regrouping, we may write Kane’s equation in the form 

N b N b 


^b;-(f-/a-wx#)J + ^fe-(E-ma)J =0 


( 2 ) 


(3) 


(4) 


k= l 


fc=i 
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The expressions in parentheses are identically zero, being Euler’s equation and Newton’s second law of 
motion for each body in the system. The vector quantities u) r and v r are the partial angular velocity and 
partial velocity , respectively. From a purely formal point of view, it doesn’t matter what uj r and v r are; 
the equation is satisfied given any bounded functions. We will discuss shortly the partial velocities in more 
detail. First, we would like to motivate that discussion by discussing the system of rigid bodies, and the 
system of equations we need to describe their motion. 

III. How Many Equations? 

A single rigid body has six degrees of freedom — three rotational and three translational. A collection 
of Nb rigid bodies may have as many as 6N b degrees of freedom, requiring up to 6 N b differential equations 
of motion to solve. For the problems that interest us, however, we connect the rigid bodies with rotary 
joints. Each joint may have one, two, or three rotational degrees of freedom and no translational degrees 
of freedom. Said another way, each rotary joint introduces three translational constraints and zero, one, or 
two rotational constraints. If N c is the total number of constraints, then there are N u = QN b — N c degrees 
of freedom remaining. 

There are several methods that we could apply to solve for the motion of the system. One method is 
to use Euler’s and Newton’s equations for each body individually, each being acted upon by the constraint 
forces (and possibly torques) arising from the joint constraints. This method requires 6 N b dynamical states 
(generalized speeds), or alternatively, N u generalized speeds and N c = 6 N b - N u constraint forces and 
torques. The joint constraints may be expressed as differential equations, 5 or as algebraic equations. 6 Either 
way, they must be solved simultaneously with the dynamics equations, yielding a total system size of 6 N b 
equations. 

Computational speed and numerical accuracy are both good reasons to look for a minimum-dimension 
system of equations. The computation time required to solve a system of equations scales as the square of 
the dimension of the system, and the forces and torques required to enforce a joint constraint may be quite 
large, driving the numerical integration timestep to be very small to keep the solution from diverging. We 
are generally not interested in the constraint forces and torques, so we seek a formulation that does not 
require solving for them. 

Kane’s method (among others) arrives at the minimum-dimension system of equations by careful choice 
of dynamical variables, and by not requiring solution of the constraint equations. By assigning as many 
independent generalized speeds as there are degrees of freedom, we partition the 6A b -dimensional space of 
system motions into two disjoint subspaces: the Ay-dimensional space of permissible motions, and the Ac- 
dimensional space of impermissible motions. All we need is something that projects Euler’s and Newton’s 
equations into that permissible subspace, and the constraints are automatically eliminated. Looking back at 
equation (4), it is precisely the partial angular velocities and partial velocities that perform that projection. 

Should it be desired to find some constraint forces or torques after all, it is possible to do so by application 
of the principle of virtual power. The reader may have been exposed to the principle of virtual work: by 
assuming an infinitesimal virtual (i.e. impermissible) motion, one may find the associated constraint equation 
required to enforce that the virtual motion is zero. The principle of virtual power is the principle that 
underlies all of Kane’s method, and it is really the same technique, using a virtual generalized speed instead 
of a virtual motion (power being the time derivative of work). By introducing, for example, three generalized 
speeds for the translational motion of a rotary joint, we may find the partial velocities and formulate equations 
of motion. Once the equations of motion are obtained, we may set the virtual generalized speeds (and their 
derivatives) to zero analytically. The virtual speeds will thus be eliminated from the equations of motion, 
restoring the minimum-dimension set of equations, and a set of three additional equations will be obtained 
expressing the desired translational constraints. We won’t discuss this technique further in this paper — we 
only mention it to emphasize the division of the motion space into permissible and impermissible (constraint) 
subspaces. The usual selection of generalized speeds spans the subspace of permissible motions, saving the 
labor of solving constraints. But we can access the constraint subspace if desired. 


3 of 13 


American Institute of Aeronautics and Astronautics 



IV. Basis Dyads, Partial Velocities and Generalized Speeds 


Consider a single rigid body. We may write the angular velocity of the body, Co, and the linear velocity 
of the body’s mass center, v, in terms of six generalized speeds: 

Co = u\bi + U 2&2 + u 3 b 3 (5) 

v = u^h 1 + ii 5 n 2 + u 6 n 3 (6) 


We may write this more compactly by introducing basis dyads*, 

M. = [fri ^2 ^3] 

K = [ni fi2 n 3 \ 

so that 

Co = Qui :3 

V = KUA:6 

Concatenating the angular and linear generalized speeds, 

Co = Qu 
v = Vu 


( 7 ) 

(8) 


0) 

( 10 ) 


( 11 ) 

(12) 


where 0= [BO] and V = [0 TV]. These dyadics are simply a compact way to express the partial angular 
velocities ancfpartial velocities. For a single rigid body, the extra notation is hardly worth it. But consider 
now a system of two (uncoupled) rigid bodies, and write expressions for their angular velocities and the 
linear velocities of their mass centers: 

uibn + W 2&12 + ^3^13 (13) 

U 4 S 21 + W 5&22 + U6^23 (14) 

Ujfl\ + U 3 fl 2 + u 9 n 3 (15) 

uioni + unn 2 + ui 2 fiz (16) 


w 1 
C 02 

Vi 

V 2 


Using the basis dyad notation, 


where in this case 


f Wi 



0 

0 

0 

1 C 02 

* 

0 

I 2 

0 

0 

1 * 


0 

0 

K 

0 

l ^2 . 


0 

0 

0 

K . 



S 

= 

li 

0 

0 

a 



z 

= 

0 0 

0 0 



g 

V 


(17) 


0 0 
0 0 


(18) 


(19) 


For any system the partial angular velocities, g, are the partials of the relevant angular velocities, Co, with 
respect to the generalized speeds, u. Likewise, the partial velocities, Y_, are the partials of the relevant linear 
velocities, v, with respect to u. For the systems under consideration, the relevant quantities are the angular 
velocities of the rigid bodies and the linear velocities of their mass centers. 

As we move from formulation to implementation, we will want to shift from vector-dyad notation to 
matrix notation, eg. from Co = Qu to to = Q.u. To do this, we will use the relation between the direction 
cosine matrix relating two frames and their basis dyads: 


B/-1A 


C A = B 1 A 


( 20 ) 


a Hughes 3 calls these vectrices. 
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Vectors and their components may also be compactly related through basis dyads. For example, let ui k — 
[uki u k 2 Wfc3] be the components of Q k in B k , and let v k = [u fc i v k2 u fc3 ] be the components of v k in N. We 
may then write 


tik = R k u k ( 21 ) 

Vk = Kv k (22) 


V. Joint Partials 


A joint imposes a set of kinematic constraints on the motion of the two bodies that it connects. For this 
discussion, we limit our consideration to rotary joints: gimballed joints with one, two, or three degrees of 
freedom, as well as spherical (ball-and-socket) joints. We need a concise and general notation to represent 
all of these. We introduce this notation by example. Consider the rotation from a body A to a body B as 
a 3-DOF gimballed joint described by a Body-3, 3-1-2 Euler rotation through angles (0 X , 0 2 , 0 3 ). We write 
the relative angular velocity of B with respect to A in terms of the generalized speeds associated with the 
joint, <7i=0i : 


cos 02 sin 03 

COS 0 3 

0 " 


f ' 

COS 02 COS 0 3 

— sin 03 

o ; 

i 


— sin 0 2 

0 



l ^ . 


(23) 


r 


a 


The appearance of signifies that the product Ter gives the components of A u B expressed in B. Analogous 
to the partial angular velocities, the matrix T relates the relative angular velocity of a joint to the joint’s 
generalized speeds. We call this the matrix of joint partials. Later in the development, we will find that we 
need the time derivative of the joint partials, f, in addition to T itself. 

Each Euler rotation sequence has a distinct joint partial matrix, with as many columns as the joint has 
degrees of freedom. For example, a Body-2, 2-3 Euler rotation through (0i, 0 2 ) has a = [gi g 2 ], and 


sin 02 

0 " 


G 2 COS 0 2 

0 ' 

cos 0 2 

0 

, . f = 

— <72 sin 02 

0 

0 

1 _ 


0 

0 


(24) 


A tabulation of w(0) relations for all Euler rotation sequences may be found in reference. 2 The joint partials 
are simply those relations in matrix form, and f is readily found by differentiation. 


VI. Two Bodies and a Joint 

We restrict our study to systems of rigid bodies connected in a tree topology, ensuring that we can choose 
as many independent generalized speeds as there are degrees of freedom. In organizing the tree structure, 
we conceptualize a building block consisting of two bodies connected by a joint. While one body may have 
multiple “children” , each joint has exactly one inner body and one outer body. This concept facilitates writing 
important kinematic and dynamic properties of the outer body ( B 0 ) in terms of the properties of the inner 
body ( Bi ) and the joint ( G k ). The angular velocity of the outer body is constructed as 

N u° = N tf i + i UJ° 

= Bui +RT k a k ( 25 ) 

The velocity of the outer body’s mass center may be written 

N v° = N v i + N uj i xr ik - N u°xr ok 

= Kvi + {Rui) x (£tn k -Kr ok ) - (RT k G k ) x {Efok) ( 26 ) 

where fi k and f ok are the vectors from the bodies’ respective mass centers to the joint. 

It is natural to express some vector quantities in terms of their components in one frame or another 
(eg. Rvi, Rxoi). It is not obvious yet which frame is most natural to express the f ’ s. We choose N in 
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anticipation of more complex systems; differences of vectors, each fixed in a different body, are expressed in 
N with equal inconvenience. 

The dynamical state vector for this two-body system is the concatenation of the inner body’s angular 
velocity components, the generalized speeds associated with the joint, and the components of the velocity 
of the mass center of the inner body expressed in N: u = [wi a\ v\] T . Writing the angular velocities and 
linear velocities for the system: 


(ft 1 

] 

1 J 

' 

f ft 1 

I = ■ 

l J 

. 




0 0 



M 

< 

i ai 

■ 

l ^ J 


N 



OJl 


o - 1 


„ Vl 


(27) 


(28) 


We may clear the basis dyads by multiplying the row by the u 2 row by and the v rows by N_ T : 


Wi 

u; 2 



U 0 0 

2 C l Tj. 0 




< 



l J 


0 

0 

u 


M 


r£ 1 w C 2 r 1 

u 





✓ 

l J 


"V* 

V 


(29) 


(30) 


where we’ve introduced /3 2 = r 2 i — rn for brevity. 

The partial velocity matrices Q, and V contain configuration-related quantities that will have to be 
evaluated at each evaluation of the equations of motion. This aside, the system velocities (angular and 
linear) are linear functions of the generalized speeds: 


{ca} = Qu (31) 

M = Vu (32) 

We present a five-body example in the appendix, to highlight the patterns in the structure of and V that 
arise. Using these patterns, partial velocities for arbitrary tree structures may be constructed by inspection. 
Carrying on, we find the outer body’s angular acceleration, 

N a° = N S i + R^iTk&k + t k a k + cj 0 x T k a k ) (33) 

and the linear acceleration of its mass center: 

N a° = x (^ x ? (i ) + V x r it - "d!” x ("i!” x f ol ) - V x (34) 

= "3* + *3* x ( w ft x r ik ) + N S i x f ik - ( w ft + BT k <r k ) x [("ft + RT k a k ) x f ok 

- N S‘ + ^(IV* + f k a k + u,x r fc o*)] x r ok (35) 


Assembling these for the system, we find that there are terms that include u and terms that don’t. We coin 
the term remainder accelerations for the latter, and denote them by a r and a r . They may be constructed 


recursively in the same fashion as a and a: 

N a° r = N 5 i r + R o (t k <T k + u 0 xr k <7 k ) (36) 

N S°r = + V x ("ff x f tt ) + x 4 - V x ("d!” x f Dt ) - ff S” r x (37) 

The remainder accelerations of the root body are zero: N a 1 r = N a} r = 0. 

The remainder accelerations allow us to write the system accelerations compactly as: 

{a} = ftu + {a: r } (38) 

{a} = Vii + {a r } (39) 
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VII. The Matrix Form of Kane’s Equation 

Equation 4 gives Kane’s equation in vector form. We have spent the intervening text developing matrix 
representations for the key factors in this equation. We have found matrix constructions for the partial 
velocities, angular accelerations and linear accelerations. The remaining factors appear in the sums over 
bodies — each body B k has a mass, m k , a central inertia matrix, I k , is acted on by active forces F k and 
torques T k , and has at any instant a gyroscopic torque term, (u x H) k . Using [m], [/], {F}, {T}, and 
{uj x H} to denote the matrices formed by assembling these for the system, we may write Kane’s equation 
in matrix form: 

« T (m - [/]{«} - X H}) + V T ({F} - [mj{a» = 0 (40) 

Substituting equations (38)-(39), and grouping on u, we obtain 

(n T [/]n + V T [m]V) u = n T ({T} - |/]K} - {u x H }) + V T ({ F } - [m]{a r }) (41) 

or, letting the brackets be implied by the context, 

(Q t IQ + V T mV) u = Q T {T ^ Ia r - u x H) + V T {F - ma r ) (42) 


VIII. Kinematic Equations of Motion 


To this point, we have focused on the dynamic equations of motion, giving scant attention to the kinematic 
relationships. For the systems of interest, the kinematics are relatively simple. As the dynamic state vector is 
arranged by the generalized speeds, u = [wi a iq] 7 ’, we arrange the kinematic state vector by the generalized 
coordinates, x = [qi 9 pi] T . 

For B i, we parameterize attitude using a quaternion, 1 q N . The rates of change of quaternion components 
are given by: 


9 l 
Q2 

1 

94 -<73 92 

93 94 — 9i 


[ : 

Q3 

1 ~ 2 

-92 9i 94 

< 

1 w 2 * 

94 , 


~q 1 —q 2 —q 3 _ 



(43) 


We have mentioned in passing the kinematic relation for a gimballed joint — representing the joint 
rotation by the Euler angles 9, the kinematic relation is simply 


9 = (j ( 44 ) 

Spherical joints are better represented using a quaternion for the joint rotation, °q i , using (43). 

The last generalized coordinate is the position of Si’s mass center in N , expressed in N. The kinematic 
relation is simply 

Pi = vi (45) 

IX. Computer Implementation 

Equation (42) is a set of N u equations in the form 


[COEF]u = {RHS} (46) 

The COEF matrix is N u x N u , and RHS is an N u x 1 array. This system is suitable for computer 
solution. The partial velocity matrices encapsulate all the coordinate transformations. Numerical solvers 
such as Cholesky decomposition exploit the symmetry of the COEF matrix, running about twice as fast as 
Gaussian elimination. 

Numerical solvers b typically require multiple evaluations of the equations of motion over each integration 
timestep. To properly assemble the COEF and RHS, the entire system configuration must be calculated 
at each EOM evaluation. Our implementation follows the algorithm presented in Table 1. 

b For our applications, a fixed-step 4th-order Runge-Kutta works fine. 
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Table 1. Algorithm for Evaluation of Equations of Motion 


For each timestep, Call 4th-order Runge-Kutta numerical integrator (RK4) 

For i=l to 4, Call Equations of Motion 

Find Joint Partials, Tfc, Tk for each joint Gk 

Map u, x to body states c Vj, ^C N , Vj, pj for each body Bj 

Find needed interbody coordinate transformations, * C k 

Find path vectors rjk and Pj 

Build partial velocity matrices 0, and V 

Find remainder accelerations, a r , a T 

Build inertia torques, {—Ia r — u x H) , and inertia forces, (— ma r ) , for each body 
Build applied torque, T, and applied forces, F, for each body, 
including active torques at joints 
Assemble [ COEF } and {RHS} 

Eliminate any locked degrees of freedom (see section X) 

Solve [COEF]u = {RHS} for u 

Find kinematic rates of change, ±=[< 71 , 9, pi] 

Return u, x to RK4 
End For 

Advance time to next step 


X. Eliminating Degrees of Freedom 

Using the generalized speeds as our dynamic variables allows us to easily eliminate a degree of freedom 
from the equations of motion. Of course, if it weren’t useful, we shouldn’t have put it in the model to begin 
with, but there are situations where it may be desirable to temporarily lock out a degree of freedom for a 
specific simulation case, or even a portion of one. For example, one axis of a joint may be locked to model 
a parking brake on a robot arm, or to model a gimbal motor failure. The form of the equations of motion 
makes it straightforward to do this without otherwise changing the spacecraft model. 

Consider the degree of freedom corresponding to Uk ■ Locking that DOF reduces the dimension of the 
subspace of permissible motions by one, consequently adding one equation into the subspace of constraints. 
We set Uk = Uk = 0, and the associated kinematic state, Xk, retains the value it had when the DOF became 
locked. The equations of motion are reduced by eliminating the fcth row and column of [COEF], and the 
ktb. row of {RHS}. 

XI. Illustrative Simulation Results 

To check our implementation, we will demonstrate a classic lesson from the dawn of spaceflight. Consider 
a spacecraft composed of a cylindrical body and four small appendages, as shown in figure 1. The spacecraft 
is spin-stabilized about its minor axis. Rotating as a solid body, linear perturbation theory shows that this 
is a stable motion. The experience of Explorer I 7 demonstrates, however, that internal energy dissipation 
can destabilize a minor-axis spin in favor of a major-axis (“flat”) spin, which has less energy for the same 
angular momentum. The internal energy dissipation was attributed to deflections of the turnstile antenna, 
excited by nutation. 

We construct a five-body dynamical model, each appendage connected to the main body by a two-axis 
gimballed joint. We introduce a spring and damping torque in each joint. To accentuate the coupling, we 
deliberately tune the joint parameters to resonate with the nutation of the system. 

The coupling is primarily in the joint axes that are perpendicular to the spin axis. When spinning about 
the major axis, only two of the joint DOFs are expected to interact strongly with nutation — we tune these 
two DOFs to resonate with the major-axis nutation frequency. The remaining six joint DOFs are tuned to 
resonate with the minor-axis nutation motion. 

Figure 2 shows the time history of the angular rates of the central body, expressed in the body frame. 
Initialized with a steady spin of 5 deg/sec about the Z axis, and a very small transverse rate (1.0E-3 deg/sec 
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Figure 1. TumbleSat Model 


about the Y axis), the nutation motion grows, eventually settling into a stable spin about the major (X) 
axis. Figure 3 shows the kinetic energy of the system, verifying that the major-axis spin is a lower-energy 
motion. The angular momentum was observed to be conserved, as expected. Not shown here, we repeated 
the simulation with all joint degrees of freedom locked, and observed that the minor-axis spin motion was 
stable under those conditions. 


XII. Conclusion 

By recasting Kane’s equations, we have shown that the partial velocities project Newton’s Second Law 
into the subspace of permissible motions, and that the partial angular velocities do the same for Euler’s 
equation. Our selection of generalized speeds is an implicit selection of which motions are permissible, and 
which aren’t. Equations of motion and constraint equations are thus related, being disjoint subspaces that 
together comprise the total 6iVb-dimensional motion space. The same tools that we use to find the equations 
of motion may also be applied to find constraint equations or to eliminate a degree of freedom as desired. 

We have used the joint-and-two-bodies building block to develop formulas that can be applied recursively 
over the system. We have used the concept of joint partials to handle all rotary joints with a common 
notation. The explicit use of basis dyads has helped us keep careful track of the relationships among the 
many reference frames involved while making the transition from vector-dyad notation to matrix notation. 

We have emphasized implementation. For systems of this complexity, it does no good to have equations of 
motion if they aren’t in a form you can solve. So we have presented not just the equations of motion but the 
algorithm for computing all the intermediate terms and assembling them for solution. The resulting system 
of differential equations is linear in the rates of change of the generalized speeds (u), and is of minimum 
dimension. 

We have limited our attention to a system of rigid bodies connected in a tree topology by rotary joints. 
This covers a broad range of spacecraft, and the concepts and notation developed here will be extended in 
future work to handle flexible bodies, translational joints, and (perhaps) more general topologies. 
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Figure 3. TumbleSat: Kinetic Energy 
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XIII. Appendix: Partial Velocity Matrices for a Five-Body Example 



Figure 4. Five-Body Dynamical Model 


To illustrate construction of the partial angular velocity and partial velocity matrices for a multi-body 
spacecraft, we consider the five-body model shown in figure 4. Using (25), we write the angular velocities of 
each body: 
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(47) 


Clearing the basis dyads, we obtain 
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(48) 


Noting that the jth row of Cl traces the path from B i to Bj, we conclude that we may write Cl for any 
tree-topology system by inspection. 

Using (26), we write the linear velocities for the mass centers of each body: 
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and, clearing basis dyads, we arrive at 


' 0 0 0 0 0 u W! 

v 2 p* N C l r 21 N C 2 Vi 0 0 0 V CTl 

' »s = A?"'? 1 0 r% 2 N C 3 T 2 0 0 U 02 

v 4 ft" c 1 o r 4 x 2 w c 3 r 2 r 4 x 3 w c 4 r 3 o u as 

Us J Lft^c 1 0 r 52 N c 3 r 2 0 r * 4 w c 6 r 4 u J CT4 

«1 

Noting the patterns in subscripts and superscripts, we conclude that we may write V by inspection for any 
tree- topology system. 
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